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Abstract 

We consider graphs made of one-dimensional wires connected at vertices, and on which 
may live a scalar potential. We are interested in a scattering situation where such a network 
is connected to infinite leads. We study the correlations of the charge in such graphs out of 
equilibrium, as well as the distribution of the currents in the wires, inside the graph. These 
quantities are related to the scattering matrix of the graph. We discuss the case where the 
graph is weakly connected to the wires. 



PACS : 03.65.Nk, 73.23.-b 

1 Introduction 

Within the field of mesoscopic physics, the interest in graphs is motivated by the fact that 
they provide simple models for networks of wires, which are most of the time sufficient to 
describe the effect of interest (such as Aharonov-Bohm oscillations of the conductance of a 
ring for example). Scattering theory plays a central role in mesoscopic physics : it provides a 
transparent formalism to study transport properties of phase coherent systems. Moreover, many 
other physical quantities can be related to scattering properties, like the current noise PJE], 
the density of states through the Friedel sum rule, mesoscopic capacitance, relaxation resistance 
[3J @] • Scattering on graphs has attracted the attention of many authors among which we can 
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Despite the scattering matrix is a global quantity characterizing the full system, some local 
information can be extracted from it. This idea has been fruitfully exploited in many works of 
Biittiker et al (see review articles [151 116| and references therein). To understand this point let 
us first consider the case of an isolated system (an isolated graph for example). In this case, the 
spectrum of the Schrodinger operator is discrete : E n , (p n (x). Let us consider a physical quantity 
described by the operator X, related to a conjugate variable / (that is X = -§jH where H is 
the Hamiltonian) . Typical examples are provided by a magnetization A4, a persistent current 



/, or the local density p(x) , which are conjugated to the magnetic field —B. a flux 1 — <p and the 
potential V(x), respectively. As it is well known, a simple way to obtain the expectation value 
of the physical quantity of interest is to compute the derivative of the eigenenergies with respect 
to / : X n = (ip n \X\(p n ) = ^jff- This result is known as the Feynman-Hellmann theorem. 

A natural question is to extend these relations to open systems that are connected to reser- 
voirs possibly in an out of equilibrium situation. As it is well known, the scattering approach 
will prove to be relevant for this purpose. The open system of interest in the present arti- 
cle will be a graph connected to some infinite wires. In this case the spectrum is continuous. 
The stationary scattering states ^^(x) describing the injection of a plane wave at contact a 
provide the convenient basis of states for the discussion. Then we can relate the quantity of 
interest (that can give a local information) to the scattering matrix X through the relation: 

\X\ tp^ ) = — (^TJj) ' ^ n ^ e case °f graphs, this idea will be made explicit in three 

cases: (i) when X — > p(x) is the local density of electrons 2 (then / — ► V(x) is the local poten- 
tial), (ii) if X — * Q is the charge of the graph (/ — > U is a potential constant inside the graph) 
and (in) if X — > is the current in an arc of the graph (/ — > 9^ is the flux along the arc). 

The purpose of our article is to study the distribution of charge and current densities in 
a graph out of equilibrium. The out of equilibrium regime is obtained by imposing different 
potentials at the external leads. A motivation for this study comes from the recent interest in 
quantum coherent devices such as Cooper pair boxes used for building charge Qubits (see [2U| 
for a review). The full spectrum of charge fluctuations is involved in the study of the dephasing 
in a Qubit perturbed by the charge fluctuations of another conductor capacitively coupled to 
the first one |15[ 1211 122 j . In the same way, current density fluctuations are a source of de- 
phasing for Qubits based on flux states. Since the formalism developed in the present paper 
provides a systematic way for evaluating the charge and current density noise fluctuations in 
a mesoscopic circuit, it might be useful for estimating the dissipation and decoherence prop- 
erties of some experimental systems of Qubits. More precisely, it was shown that transition 
rates of a two level system weakly coupled to a quantum environment are directly related to 
the unsymmetrized correlator : see for example |23| where the roles of the negative and the 
positive part of the spectrum of the unsymmetrized correlator are studied. On the other hand 
the relaxation and decoherence rates are related to a symmetrized correlator [SO]. Correlators 
are more directly accessible in noise measurements : in a recent work, Gavish et al. |24| I25| 
proposed a description of the full measurement chain for the current noise of a mesoscopic 
sample. In this work, the unsymmetrized correlator is involved in excess noise measurement. 
Finally one should mention that experimentalists are now able, using photon assisted tunneling 
in a superconductor-insulator-superconductor tunnel junction, to measure the unsymmetrized 
current noise correlator in quantum mesoscopic devices j^H]- The question of which correlator 
(symmetrized or not) to consider depends on the question of interest. Therefore we will consider 
in the following the unsymmetrized correlator as the fundamental object 3 . 

In this paper, electron-electron interactions will not be taken into account. However, even 
if this limits the applicability of our results, we recall that they can be taken into account in 
a mean field Hartree approximation within the scattering approach, as it has been developed 
in several papers by Buttiker and collaborators [H| H| l2"7] (see also |28j for a review). In this 
framework the charge (or the current) contains two contributions : a bare contribution (injected 

The variable conjugate to a flux line threading a loop of a planar graph is the current flowing through the 
semi infinite line issuing from the flux [HI (see also [T5|V Th is can also be easily understood in the 2-dimensional 
plane |19| . 

2 ln this case {'tp^ \p( x )\ Ve ) = ve (x)^^' '( x ) ls an °^ diagonal element of the local DoS |15l II3| 
3 There are indeed two unsymmetrized correlators depending on the order chosen. But for a system in a 
stationary regime, they can be simply related as shown in appendix 1X1 
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charge) and a contribution from screening. Screening affects ac transport or finite frequency 
noise. It is not the purpose of our article to consider such interaction effects; in other terms we 
will focus on the bare contribution of the charge and its relation to the scattering matrix of a 
graph. 

This paper is organized as follows: first of all, the basic formalism lij necessary for the 
discussion is recalled. Then the charge distribution inside the graph is analyzed in details. 
The first and second moments of the total charge are related to the scattering matrix. Finally 
analytic expressions for the full spectrum of charge fluctuations are provided. 

In a second part we will show how currents inside the graph can be related to the scattering 
matrix. It was shown in |17j that the persistent current can be related to the derivative of 
the Friedel phase with respect to the magnetic flux. In this work the possible generalization 
to an out of equilibrium situation was not considered because the authors did not identify the 
different contributions of the various scattering states associated to the different leads. These 
contributions were identified later in |29j . Moreover, in Taniguchi's work, a formula relating the 
current-current correlations and the scattering matrix was proposed. Despite the contributions 
of the scattering states to the correlator were given, it is still not sufficient to study current- 
current correlations in an out-of-equilibrium situation (this point will be made clear later). Our 
results go beyond this limitation and provide the generalization of Taniguchi's result. 



2 Basic formalism : scattering matrix 

We consider the Schrodinger operator — D x + V[x) on a graph, where D x = d x — iA(x) is the 
covariant derivative (we choose units : h = 2m = e = 1). The graph is made of B bonds (af3), 
each being identified with an interval [0, l a p\ G M. We call x a p the coordinate that measures 
the distance from the vertex a. The Schrodinger operator acts on a scalar function (p(x) which 
is described by B components f( a (3)(x a /3), one for each bond. The bonds are connected at V 
vertices. The adjacency matrix a a p encodes the structure of the graph : a a p = I if (a/3) is a 
bond and a a p = otherwise. 



Vertex formulation 

Let us first assume that the wave function is continuous at each vertex. This allows to introduce 
vertex variables ; we denote (p a = if (a) the function at the vertex a. The continuity condition 
reads : ip( a p)( x ap = 0) = y? a for all vertices /3 neighbours of a. A second condition is added 
to ensure current conservation at the vertices : Y^B a aP^x t -P(a,p)( x aP = 0) = X a (p a where the 
sum over (3 runs over all neighbouring vertices of a due to the presence of the adjacency matrix. 
\ a is a real parameter. The requirement of continuity of the wave function imposes a special 
scattering at the vertices : in particular, the transmission amplitudes of a plane wave of energy 
E = k 2 between two leads issuing from the same vertex of coordinence m a = Y^p a af3 ar e all 
equal to 2/(m a + i\ a /k). 

The wave function on the bond (a/3) is : 

<P(»0)M = e ix ^/ 1 ^ (<p a f a p{x aP ) + ipp e- i6 <*ef Pa (x a pj) (I) 

where a p is the magnetic flux along the bond (a/3) (the vector potential is A a p = 9 a p/l a p). The 
two real functions f a p{x), f/3 a (x) are the two linearly independent solutions of the Schrodinger 
equation [E + d x — Vr a p\(x)]f(x) = on the bond satisfying boundary conditions : / Qj a(0) = 1, 
fapilap) = 0, f/3 a (0) = and f(3a{laf}) = 1- These two functions encode the information 
about the potential on the bond. For example, in the absence of potential, V(x) = 0, we have 

. , x _ &\xik{l a(1 -x afi ) 
Ja/3\x a /3) — s inkl a/3 
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The graph is connected to L leads. Each lead is a semi- infinite line plugged at a vertex of 
the graph, with a coupling parameter w a £ R (see figure^). The introduction of these couplings 
allows to go continuously from an isolated graph to a connected one. The precise physical 
meaning of these parameters is given in ^U]. In particular, the transmission amplitude through 
the box between the graph and the lead is 2w a /(l + w^). We introduce the L x V matrix W : 

W a/3 = w a 5 a /3 (2) 

where a belongs to the set of vertices connected to leads and (3 to the set of all vertices of 
the graph. This matrix encodes the information about the way the graph is connected to the 
external leads. 

The scattering matrix S is a L x L matrix describing how a plane wave of energy E entering 
from a lead is scattered into the other leads by the graph. It is given by 



where the matrix M is : 



M aP (-E) = — \5 a p 



a ^T—( a ) 



+ a a ^(P)e iS A . (4) 



dx a(3 



Note that for E > 0, this matrix is antihermitian : = —M. It can also be related ^Uj to 
reflexion and transmission coefficients describing the potential on each bond : 

M a0 (-E) = 6 a 8 (i^ + y a a „ j 1 ~ ^jS 1 + r » a \ + W ^ a 1 
~^ ' ~" \VE ^ ^(l + r^)(l + 

~ a cxfij\ i \f\ i \ — 7 — 7 — • y b 

(1 + r a/ 3)(l +rp a ) - t a pt Pa 

The expressions ()4I">|) , together with @ , generalize results known in the absence of the potential 
0. 




Figure 1: Example of graph. The boxes represent the couplings between the infinite leads and 
the graph. 

The vertex formulation that we have just recalled is rather efficient mainly because vertex 
matrices are rather compact. However, as we have noticed, it corresponds to a particular choice 
of vertex scattering which does not describe all allowed relevant physical situations. In the most 
general situation it is not anymore possible to introduce vertex variables and one has to use the 
arc formulation that will now be briefly described. 
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Arc formulation 



An arc is an oriented bond. On each arc i we introduce an amplitude Ai arriving at the vertex 
from which % issues and an amplitude Bi departing from it (see figure [2J) • Equivalently, the 
wave function ipi( x ) on t ne bond is matched with Aie~ lkx + Bie lkx at the extremity of the arc. 
It is clear than we have to introduce L such couples of amplitudes, one for each external lead. 
These external amplitudes are gathered in L-column vectors A ext and B cxt . By definition the 
scattering matrix relates these amplitudes : B cxt = T,A ext . On the other hand we must introduce 
two couples of amplitudes Ai, Bi per bond of the graph, i.e. one couple per arc. We gather 
these 2B amplitudes into the column vectors A mt and B mt . Finally we group all amplitudes, 
internal and external, in two 2B + L column vectors A and B. 

The scattering by the bonds is described by a matrix R coupling reversed internal arcs : 
A mt = RB mt . The matrix element between arc i and j is given by : 

Rij = ri5ij + tjSij (6) 

where i designate the reversed arc. and ti are the reflexion and transmission coefficients 
describing the scattering of a plane wave by the potential of the bond (i). The scattering at the 
vertices is described by a matrix Q coupling arcs issuing from the same vertex : B = QA. If the 
basis of arcs is organized as {internal arcs, external arcs}, the matrix Q can be separated into 
blocks: 



Q 

The scattering matrix reads : 



Q' mt 


Q T ' 


Q 


Q ext 4 



(7) 



£ = Q ext + Q (ijt - Q int ) _1 Q T . (8) 

For more details, see |lflj . Historically, the arc approach has been followed in many works, like 
|3U1 1311 132j since it is the most natural approach. It has been formalized more systematically in 
|33| without potential and in ^U] in the most general case. 

How can we express the wave function inside the graph within the arc formulation ? In 
this case, the appropriate basis of solutions of the Schrodinger equation on the bond [E + d^ — 
Vr a p\(x)]f(x) = is not anymore the functions f a p{x) and fp a (x) introduced above, but the 
couple of stationary scattering states 4> a p(x) and 4>f3 a (x) associated to the potential Vt a p\(x) on 
the bond (a/3). If we imagine that the potential Vr a p\(x) is embedded in R, then the function 
<j> a 3(x) is the scattering state incoming on the potential from the vertex a and is matched out 
of the bond to : <j>ap(x) = e lkx + r af3 e~ lkx for x < and 4> a p(x) = t al 3e lk( - x ~ la ^ for x>l a p jUJ]. 



Figure 2: The internal amplitudes associated to the arcs a/3 and /3a. 
Then the component of the wave function (p(x) on the bond (a/3) reads : 

<P{ a p)( x otp) = B a/ 3 4> a p(x aj3 ) + Bf3 a <f)p a (x a p) . (9) 
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3 Charge of the graph 



This section is devoted to the study of the charge distribution in the graph. Our discussion will 
focus on the average and correlation of the total charge of the graph. The average charge and 
the zero frequency charge noise can be related to the graph's scattering matrix. These relations 
provide an extension of the Feynman-Hellman theorem for open systems in an out of equilibrium 
situation. Then, we present a detailed study of the charge noise at finite frequency, emphasizing 
the effect of the non equilibrium regime. For simplicity, we shall work at vanishing temperature. 

It is convenient to use the language of "second-quantization" and introduce the field opera- 
tor : 

L 



a=l 



#M) = V / dE^(x)a a (E)e-^ t , (10) 



where a a (E) is the annihilation operator associated to the stationary scattering state ^^(x) 



(a) 
E 



corresponding to a plane wave of energy E injected from the lead a. Note that ip{x, 0)| 
ip^ (x) | vacuum ) . 

Studying the charge distribution for graphs with localized states jHU E3 would re- 
quire taking into account the contribution of the discrete spectrum in the field operator : 

Sf=i Pnjix) a-n,j e~ lEnt (the function tp n j(x) is an eigenstate of energy E n localized in the 
graph and thus normalized to unity in the graph, and j denotes a degeneracy label). However, 
such a situation in not generic but arises from symmetries of the graph. For this reason, we 
shall not discuss it here. 

In a non equilibrium situation, the quantum statistical average gives : {a\ x {E)ap{E')) = 
5 a /35(E — E') f a (E) where f a (E) is the Fermi-Dirac distribution function giving the occupation 
of the scattering states coming from the lead a. 

Charge operator. The charge operator is : 

Q(t)= dx^(x,t)ip(x,t) . (11) 

J Graph 

We introduce its matrix elements on the shell of energy E : 

p^\E) = {^\Q{t)\i>f). (12) 

Since the spectrum is continuous, these matrix elements have the dimension of a density of states 
(DoS). They can be related to the scattering matrix : 

^) (B) ._J_( E t|^, (13) 

where U is a constant potential added inside the graph only (the variable conjugate to the charge 
of the graph). Instead of differentiating with respect to some additional background potential 
U, it is also possible to relate it to the derivative with respect to the energy : 

n ~ L> dl *' (l) " if{x) = h + > " E, ')„, • < 14 > 

Note that X)aP K-^) ^ s ^ ne of the graph, i.e. the local DoS integrated inside the graph. 
These relations are proven in appendix O 

Average charge. The average charge : 

=E [°° dE UE) I dx\^\x)\ 2 = Y, r dEf a {E)p^{E) (15) 

JO ./Graph J 
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involves the injectivities p(x,a;E) = \ip^\x)\ 2 , which are the contributions to the local density 
of states (LDoS) coming from the scattering state tp^\x). The contribution p^ a,a \E) of the 
scattering state to the DoS of the graph is weighted by the occupation Fermi factor in 
the lead a. This illustrates the necessity of the concept of injectivities, emissivities, etc, in the 
context of out of equilibrium systems EE] (see also JH] for a discussion in the context of 
graphs). 

Charge correlation function. The charge correlation function is defined as : 



+00 



S QQ (u) = d(t - (<Q(t) Q(t')) - (Q(t))(Q(t'))) e-(*-*') . (16) 
J —00 

In appendix the relation between this unsymmetrized correlator and the other one (Q(t) 
and Q(t') in reverse order) is clarified. Definition IJ16JI matches with the one used reference [22] 
which contains a detailed discussion of the relation between the unsymmetrized correlator and 
transition rates in a two level system linearily coupled to the charge operator. 
Using the relation (atapala,,) — (a^ap)^ 

O'v) — "03/ 8/3n /a(l - fp) we obtain : 



(17) 



S QQ (u)=2nJ2 I™ AEUm-UlE + u)] [ dx ^ (x)* ^Jx) 

a a J J Graph 

Only the zero frequency correlations involve the p^ a ^ (E) : 

POO 

S QQ (u = 0)=2nJ2 / dEf a (E)[l-ME)]p^(E) P ^ a) (E). (18) 

a,P J ° 

Using equation (|14jl the zero frequency noise of the total charge is related to the scattering 
matrix of the graph. 

Charge fluctuations. The charge fluctuations at a given time involve the integral of the full 
spectrum : 

Q2 = (Q(t) 2 ) - (Q(t)) 2 = ^J^ S qq (lo) . (19) 
In terms of the stationary scattering states, we get : 



12 = 

a,p 



/*oo p 2 

&EdE> UE)[l-f p {E>)} / dx^\xT^}{x) 

„, a JO J Graph 



(20) 



Weakly connected graphs 



To go further let us focus on the case of graphs weakly coupled to the leads {w a — > 0) . Note that 
we do not consider charging effect in the following (Coulomb blockade) which is important if the 
capacitance describing the Coulomb interaction between the leads and the graph is small (see 36, 
for a review article). A description of such effects would require a different approach. However, 
in the neighbourhood of the Coulomb peak, a description within the scattering approach can 
be sufficient to describe transport, like it has been done very recently in [HZj to analyze Fano 
profile measurements in a ring with a dot embedded in one of its arm. 

If w a — > the decomposition of the scattering states over the resonances (levels of the 
isolated graph), derived in appendix El can be used 
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Here <p n (x) denotes the wave function of the eigenstate of energy E n of the isolated graph, 
normalized to unity in the graph. From this expression we get : 



where F nia = y/E^w^\ip n (a)\ 2 is the contribution of the contact a to the resonance width 
r„ = £ a r n , a . The phase is given by e 1 ** = gfgg^§ - 

Average charge 
Equation (fT5|) gives : 



(Q(t)) * x 



f - /.<*> E - E E %f (i -an ^ + I) 



Of " " n K • na 



(23) 

where the sum over n runs over the energies of the resonances (energies of the isolated graph) . V a 
is the potential at contact a. This equation was derived in jSS] by tracing out the lead's degrees 
of freedom. Since the average charge is the sum of contributions of the various levels, we can 
consider only one level E n . If the level is below the potentials, E n < Vr < Vl, the occupation of 
the level is 1. On the other hand, if the level E n is between the potentials, Vr < E n < Vl, and 
far enough from them (on the scale T n ), it gives a contribution r n> L/T n to the average charge, 
which simply expresses that only the left scattering state is contributing to the occupation of 
the resonant level. 



Charge noise at finite frequency 



Let us now discuss the finite frequency structure of the charge noise for weakly connected graphs. 
Equation (|TT)) requires evaluating 



Graph 



TT 



r E n w 

(«K(/3) 



TT (E E n iT n )(E + u-E n + ir n ) 



(24) 



Let us keep only the diagonal elements in the double sum. This diagonal approximation is valid 
in the limit of narrow resonances (T n <C \E n+ i — E n \) since the energies E and E + co are then 
compelled to be both in the neighbourhood of the level E n . Then the correlation appears as a 
sum of contributions of the different energy levels : 



n 

The contribution of the level E n reads : 

sg,M = ^Ej dE um - uiE+»)\ Tn ' a/7T 



(25) 



r n ,/3/v 



Performing the integrals leads to 

1 



(E - E n y +Tl(E + uj- E n y + vi 



(26) 



(27) 
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where 6{u) is the Heaviside function and 

V n — E r , 



A(V a , Vp;u) 



arctan 



+ arctan 



arctan 



Vp — E n 



V a +u - E n 



arctan 



Vp — u — E n 



+ Vn ln \{V a +u- E n f + T 2 n ][(Vp -uj- E n f + Tl 

UJ 



[{y a -E n Y + Tl][{Vp-E n f + Tl\ ■ 
For the particular case of a two terminal geometry the noise reads 

Sqq(<») = J r3l + J /4r2 [rl L e(u)A(V L ,V L ;uj)+Tl R e(uj)A(V R ,V R ;uj) 
+T n ,iXn,R B(uj + V) A(V L , V R - uj) + T njR r niL 6(uj - V) A(V R , V l ;uj)} , 

where V = Vl — V R > is the voltage drop. 
• Equilibrium case : Vl = V R = 0. 



Vr. 



t-'n 



Vr 



wire L Graph wire R 



(28) 



(29) 



2TrT n 1 + w 2 /4r2 



u - E n \ ( ' uj + E r , 
arctan — — + arctan 



IV \{uj-E n f + Tl][{uj + E n f + Tl] 
[El + Tlf 



(30) 



We consider the case of narrow resonances where T n is the smallest energy scale, then if the 
frequency is smaller than \E n \, the contribution is zero, but if u is sufficiently large to excite an 
energy level u > \E n \, we get a contribution : 

Mr. a ~ 1 g H 

2r n l + u; 2 /4r2 




W«)*^T7^**< , . (31) 



Obviously, the transition between the two results is not sharp but occurs on a scale T n . Practi- 
cally all the noise power is concentrated at low frequencies as shown on figure El 

Note that this contribution is independent of the fact that the level is occupied (E n < 0) 
or empty (E n > 0) since it is an even function of E n . At Vl = V R , the low frequency charge 
noise can be understood using a classical stochastic model describing the relaxation process of 
an electron (or a hole) with lifetime l/2T n . 

• Non equilibrium regime : Vl ^ V R . 

Zero frequency limit. Then only the term with A(Vl, V r ; 0) contributes : 

Sqq(lo = 0) ~ — n,I l „ n ' L I arctan ( — ^— — -J + arctan ' — — — 



Tn(VL-E n ) T n (E n -V R ) 
{V L -E n y + Tl + (E n -V R r + T2j ^ '■'-> 
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Figure 3: Equilibrium noise in terms of —10 < oj/Y n < 10 and of the energy of the level 
(—E n )/T n (this is equivalent to vary the energy or the chemical potential). 



Each energy level brings a contribution only if it is between the potentials (Vr < E n < V£) and 
far enough from them (on the scale T n ). In this case, the zero frequency charge noise is given 
by: 

= 0) * . (33) 

n 

We recognize the factor r nj ^r n) £ characteristic of partition noise : if one of the couplings vanishes 
(Xn.B, = or F nj L = 0) the occupation of the level is either or 1 and does not fluctuate. These 
fluctuations are a signature of the non equilibrium regime of the mesoscopic circuit. The charge 
fluctuations at fixed time can be obtained by integrating the noise spectrum (j!9j) . Using the 
approximate expression 1)21)1 in 1)20)1 we obtain 



X — * X — v Fn aF n a I 1 V a — E n 1 \ / 1 E n — Vff 1 \ 

52 * \ L -r^- arctan + 2 ) arctan + 2 J (34) 

In the two leads case, only the levels between the two potentials bring the contribution : 



l\i IAn ,R /or\ 
Q2 ^ 2 • (35) 



Note that this result cannot be simply infered from the current shot noise. Near a resonance, 
the transmission probability through the graph is T(E) ~ j^-^^j^- . The average current in 

the lead is given by the Landauer formula (J) = ^ fy^ dET(E) whereas the current and the 

shot noise by S n {uj = 0) = ± dET(E)(l - T(E)) [TJ H ESI . If only one level E n lies 

r r 

between the two potentials, we obtain in this non linear regime HO] : (J) ~ 2 n 'f n ' L and 
S H (u J = 0)^2 T -^^(Tl iR + Tl iL ). 

Finite frequency noise. Let us choose the origin of the energies in such a way that : Vr = 0, 
Vl = V > 0. Four energy scales must be considered : E n , T n , V and u>. Several regimes can 
be observed according to the frequency. To help the discussion we neglect the smallest scale, 
supposed to be T n , as we did above. 

(i) Let us first discuss the case of a fully occupied level : E n < Vr = 0. 
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Graph 



At small frequencies, correlations are roughly zero. When uj reaches VR — E n , contributions from 
the second and third terms appear, while all terms contribute for uj larger than Vl — E n . To 
summarize these three regimes : 

An), \ _J 1 

where we have factorized the equilibrium result. 

In the second regime Vr — E n < u < Vl — E n , the energy to is sufficient to excite the state 
originating from the right reservoir but not from the left reservoir. This is the origin of the 
ratio T n pi/T n . In the third regime, both reservoirs contribute to the noise. At fixed non zero 
bias voltage, this leads to a double peak structure in terms of the uj which corresponds to the 
threshold for creating electron-holes pairs involving the left and right leads (see figure . At 
small V these two peaks tend to merge into a single more pronounced one. At large V the second 
peak occurs at a larger frequency and is less pronounced because of the Lorentzian factor. 



( if uj<V R -E n 

^ if V R - E n < to < V L - E n 



(36) 



1 



if V L -E n <co 




Figure 4: Non equilibrium noise in terms of —20 < Lo/T n < 100 and of the voltage drop V/T % 
for a fully populated level (Vr = E n + 50r n ). 

(ii) Let us now discuss the case of an empty level : E n > Vl- 



V L z 



L -'n 



Vn 



Graph 



The physical picture can be obtained using a hole picture. Three different regimes can also be 
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distinguished : 



4£h 



2T n 1 + w 2 /4r2 



x < 



r o 

r„ 

1 



if u<E n -V L 

if - < W < " Vr 

if - Vr < W . 



(m) Finally we consider the case of a level between the two potentials Vr < E n < Vl- 



Vlz 



m 



Vr 



Graph 

Considering the case where the level is closer to Vl than to Vr, we obtain 



2r n i + w 2/4r2 



X < 





— T 1 

T 2 



T 2 " 



1 



if U <-E n + V R 

if -E n + V R <to< -V L + E n 

if -V L + E n < u < V L - E n 

if Vl - E n < to < E n - V R 

if E n -V R <u 



(37) 



(38) 



The fluctuation spectrum is symmetric in the interval centered around uj = with width of 
order V . The main contribution to the noise appears at low frequency as can be seen from figure 
121 The multiple plot (figure |HJ) shows how the low frequency peak develops when Vl crosses the 
energy level. 




Figure 5: Non equilibrium noise in terms of —40 < uj/T n < 40 and of V/T n . The right lead 
chemical potential is fixed to — 10r n . Case (in) corresponds to V > WT n and exhibits an 
important low frequency noise whereas V < 10T n corresponds to case (ii). 

Interestingly, correlations are proportional to the partition factor T n ^ R T n L only for small 
frequencies \oj\ <C V (or large time scales t 3> ^-/V). For large frequencies |cj| 3> V, the partition 
factor does not appear. The high frequency part of the charge fluctuation spectrum is insensitive 
to the fact that the system is out of equilibrium. In this limit, the equilibrium result (|3*Tj) is 
recovered. 
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-10 







10 



co/r 



Figure 6: Multiple plots of non equilibrium charge noise in terms of —20 < u:/T n < 20 for 
different values of V/T n . The right lead chemical potential is fixed as on figure (a) V/T n = 4, 
(b) V/T n = 8, (c) V/T n = 9 and (d) V/T n = 14. 

4 Currents inside the graph 

A possible way to probe a mesoscopic device is to attach some leads to it, through which some 
currents are injected. Some information can be extracted from transport or noise properties : 
average values and correlations of currents in the external leads. All these properties can be 
related to the scattering matrix (see [40] for a review). If one is now interested in local infor- 
mation on the system, like the measurement of a persistent current, a natural way would be to 
introduce some local probe. However, as we have recalled in the introduction, local information 
can also be extracted from scattering properties. Here, we investigate the currents in the inter- 
nal wires, and show the relation to the scattering properties. The starting point, exposed in the 
introduction, is the relation between the current in a wire and the derivative of the scattering 
matrix with respect to its conjugate variable, the flux in the wire. This idea comes from |17j and 
has been elaborated further in |29| to include derivation of correlations of the current density. 
Here, we focus on the case of graphs in the context of which we will generalize these previous 
results to the non equilibrium situation. 

4.1 Current in a closed graph 

First we derive an expression for the current density in a closed graph. It is convenient to 
introduce the spectral determinant of the Schrodinger operator 

5( 7 ) = det(-A + V(x) + 7 ) = \{{E n + 7 ) (39) 

n 

where 7 is a spectral parameter. The set of E^s is the spectrum of the Schrodinger operator 
on the graph. It was shown in |41[ 1421 I43| I44j that the spectral determinant, which is the 
determinant of an unbounded operator, can be related to the determinant of a finite size matrix. 
In the vertex approach, the formula (|100j) involves aVxV matrix, whereas in the arc language 
the spectral determinant involves a 2B x 2B matrix : S(—E) oc det(l — QR) |42l 144] . Note 
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that for a closed graph, the vertex scattering matrix Q has the same dimension as the bond 
scattering matrix. We introduce the current density j(E) associated to the states in the interval 
[E, E + dE[. The current density in the arc a is : 

F)W 1 f) 

j a (E) = -Y^5(E- E n )-^ = - Im — In S(-E + 10+) , (40) 
where a is the magnetic flux along this arc. 

Example : Consider a closed ring of perimeter / threaded by a flux 8. Its spectral determinant 
is 5(7) = ch( v /7 1) - cos(0) [HI. We write : 7 = -A; 2 + i0+, then ch^Z) = cos(/d) + i0+ sin(/cZ) 
and we get for the current density in the ring : 

j(E) = -sm9 sign(sin kl) 5 (cos kl - cos 9) = ^5{E - E n )I n , (41) 

n 

with I n = —dgE n where E n = (2nir — 6) 2 /l 2 . 



4.2 Current in open graphs 

Now we consider a graph connected to infinite leads. 
The current operator is 



j(x,t) = j ip j (x,t)D x 4i(x,t) -B* x ft(x,t)ip(x,t) (42) 



where D x = d x — \A(x) is the covariant derivative. 
We introduce the current matrix elements 

jfc® (E) = ( \J(x,t)\ j>f > for x G iiv (43) 

which can be shown to be independent of the coordinate x along the arc (only if the two states 
have the same energy). This matrix element can then be computed at the vertex \i [x = 0) : 

j^XE) = \ (^D^M -D^gjfci)^) . (44) 
The quantum and statistical average of the current operator in the arc \iv gives : 

<V = (J(x e fiu,t)) = W dEf a (E)jft a \E) . (45) 

The correlations (unsymmetrized in time or frequency) between the currents in the arcs \iv and 
fi'u' , defined as 



= / d (*-*') ((j(x,t)j(x',t>)} - (J(x,t))(J(x',t'))) e^-t) , (46) 
J— 00 

for x £ [iv and x' G /uV. They can be rewritten at zero frequency as : 

Sj^Au = 0) = 2vr^| dEf a (E)[l - f p {E)\ jjgA \E) jffi »(JB) . (47) 
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Vertex formulation 

Now we look for a relation between the current density and the scattering matrix. We start 
from the expression of the scattering state in the arc jiu : 

^V) (*) = e^/W (W Uu(x) + ^ e- ie -/^(x)) . (48) 



Then : 

= 1 ^ {v) e -i^ $p + g^) e iV ^ (49) 

where we have used the fact that the Wronskian of f^v(x) and f v ^{x) reads : ^r^(p) = ~^r^( L ')- 
From the definition of the matrix M, we see that : 

Only the two elements M^ u and M Ufl depend on the flux 9^, then : 

where VP is the V x L-matrix that gathers the values of the L stationary states at the V vertices : 
*Aia = ffi. This matrix is [TOUEl : 

* = jswrwm wT ■ (52) 

We can rewrite the current density in terms of matrices M and W : 

= -1 (w I ™- I wA (53) 

J ^ iTT \ —M + W T W d9^ M + W T W J a/3 { } 

Our aim is now to find the relation of this expression with the scattering matrix. We use the 
relation ^A(r])~^ = — A(rj)~ l d A(rj)^ 1 that gives the derivative of the inverse of a square 
matrix A{rf) depending on a parameter rj. In the expression Q only M depends on the fluxes, 
it follows that : 

^ = -2W l - ™- l - (54) 

dO^ M + W T W dO^M + W T W ' { } 

If we multiply this expression by £t from the left, it replaces the M in the left fraction by —M. 
We conclude that the off-diagonal elements of the current density reads : 

- 1 ■ w> 

Note that this result is reminiscent to the one obtained by Taniguchi in j^S] who derived some 
relation between the scattering matrix and the "current density", i.e. the diagonal elements 
(a = /3) of (jSSJ). 
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Arc formulation 



Let us now reformulate the previous demonstration in the arc language which allows to consider 
the most general case. 

As explained in the introduction, the wave function on the bond (nv) can be expressed as : 

VVO^M") = B^ v ^(x/ju,) + B Vfl (punixnu) (56) 

where 4>fj, u {x^u) and fyv^ix^y) are the left and right stationary scattering states for the bond 
potential Vn u ,\(Xf jU/ '). Using the expressions of these functions at the extremities of the bond 
given in the introduction, we get the derivative of the wave function at the vertex fi : 

DxVVi/) (a*) = ik [B^ u (l - r M „) - B Vfl t Vfl ] . (57) 

We call Bj%) the internal amplitude corresponding to the stationary scattering state ^^\x). 
These amplitudes are obtained by solving the equations B = QA and A mt = RB mt with external 
amplitudes A cxt describing the injection of a plane wave on the lead arriving at vertex a : its 
components are A^ exi = .) , & arca where "area" designates the arc related to the lead issuing 



from a. Then the amplitude 

^1 - Q int J R)- 1 Q T l (58) 



is related to the matrix element between the internal arc fiv and the external arc "arc a" . 
After a little bit of algebra, we get for the current density on the arc [iv : 

Aa,/3) _ ou |2 ft(fi) i 6(a)*+* r £}(£) _ f>{a)* f * r 5(/3) _ 5(a)*|+ PrO 3 )! (^Q) 

We have used t*^ v r vtl = —r^ u t Ufl coming from the unitarity of R. In the bond scattering matrix 
R, only the transmissions depend on the magnetic fluxes : t^ u oc e l6,M ". It follows that, in the 
matrix 4^—R, only the 2x2 block related to the arcs fiu and v\i is different from zero. It is 
given by : 

Then it is straightforward to see that : 

3^ = -2^^* (^R) Bf (61) 



1,3 



where the sum over i,j runs over the 2B internal arcs. Using the expression l|58jl for the 
amplitudes, we obtain : 



J 



fcfi = -±(q*{1-R^Q^)-^R(1-Q^R)^qA (62) 
2ivr V de^ J a fj 

= -J- ( EtQ(i?t - Qint)-l^ (jR t _ Q^-IqA (63) 
2ivr V dd^ J a j3 

- U^&d-^h-U^h- (64) 

We have recovered the formula (|55|) within the arc language 4 . This demonstrates that equation 
(f53|) applies to the most general situation, as expected. 

4 We have used the relation E t Q(-R t — Q lnt ) _1 = Q*(l — R t Q lnt coming from the unitarity of the scattering 
matrices 
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Average current and current correlations in terms of the scattering matrix 

The average current can we written, as could have been guessed from the general discussion of 
the introduction : 

J ~ = sf (* at) ' (65) 

The correlations of currents at zero frequency rewrite in terms of scattering matrix : 

w - o) - 4g - « (* ^ • <«> 

If jj,u = /iV this gives the noise of the persistent current. At equilibrium, all potentials are equal 
f a {E) = f(E) Va, and we recover an expression reminiscent of the one given in |29| : 



Sj^Aco = 0) = ±-J dEf(E)[l - f(E)} Tr j 



dS d£t 



(67) 



In his work, Taniguchi identifies the contribution ( , dSt ) of a given scattering state to 

this trace. However we see that it is not sufficient to go back to the expression (|66|) describing 
the non equilibrium situation. 

4.3 Gauge invar iance 

Since many formulae involve the fluxes along the wires, it is important to discuss how a gauge 
transformation would affect them and to check that all the measurable quantities are indeed 
gauge invariant. A gauge transformation changes the vector potential according to A{x) — > 
A(x) + d x x(x), where x( x ) is a scalar function. The magnetic flux 6^ u along the arc fiu is then 
modified according to : 

— *■ ®'nv — Ofiv + Xm - Xu (68) 

where \n = x(m) is the value taken by the function at the vertex fj,. In the vertex approach, we 
immediatly see from its defintion that the matrix M is changed as : 

— M'^ u = Mfxv e ix *~ ix " . (69) 

We can write M' = lAMlA^ where the diagonal unitary matrix reads : lA a p = 5 a pe lXa . Since 
W T W is also diagonal it is clear that (±M' + W T W)- 1 = U(±M + W T Wy 1 W . The scattering 
matrix changes in the same way : 

S Q/3 — K/3 = e**--*" . (70) 

From H53|) or (|55jl we see that the matrix elements of the current operator pick up a phase 
through a gauge transformation 



J 



Nevertheless, the average current (|45j) and the correlations ()47p are gauge invariant, as they 
should. 
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4.4 Weakly connected graphs 

As we did for the charge distribution, it is interesting to consider the case of graphs weakly 
connected to the leads (w a — ► 0) for which interesting results can be derived. The starting 
point is again the expression 1)21(1 of the scattering state near a resonance (when E is close to 
an eigenenergy of the isolated graph). Using this relation, the current density matrix element 
can be expressed as 

MP)(E,E') =(^\j( x efiu,t = 0)\^)) (72) 



>liv V-'i J -' ) — WE \" v" ^ P"i " ~~ YE' 

fRn W a ip n (a) WpiPniP) 



(73) 



E,E'~E n vr E - E n - \T n E> - E n + iY n 
where J" is the current in the arc \xv associated to the eigenstate ip n {x) of the isolated graph : 

= ~ il Pn( x ) VxVnix) + cc. for xe fiu. (74) 

Note that in principle the current matrix element (|72jl between scattering states of different 
energies depends on the coordinate x, however in the weak coupling limit, since the two scattering 
states are proportional to the same eigenstate ip n (x), the matrix element becomes x independent. 

Equation (|72|) shows that the calculation of the average current is very similar to the calcu- 
lation of the charge (|23|) : 



(75) 



The contribution of the resonant level can be written : 

Jt^iLiQit)}^ (76) 



where (• • -p n ' designates the contribution of the resonant level n. 
Similarly we obtain for the correlations : 



For example, for a situation with two contacts with a potential drop V and only one resonant 
level contributing we get : 



and 



■V * (78) 



V,v(0)^^v^ (79) 
4.5 Example 

Let us focus on the simple example of a ring with two leads (see figured). 
The scattering matrix of the ring reads : 

2 / \w\sm.kl + w\wls a s b iwiw 2 {s b e~ l9a + s a e ldb ) \ 

S = "I + -H n (80) 

S y iu>2Wi (sj,e a + s a e b ) \w\ sin kl + wfwls a Sb J 

where 9 a and 9 b are the fluxes of the two arcs and 9 = 9 a + 9 b the total flux threadening the 
ring. We have denoted s a)b = sinA;/ aj fe. 

S = Sa s b det(M + W T W) = 2(cos 9 - cos kl) + i(w{ + u|) sin A;/ + w|s a s 6 (81) 
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Figure 7: A ring with two arms of lengths l a and lb, threaded by a flux 9 and coupled with 
two leads, with coupling parameters W\ t 2- The boxes represent the tunable couplings, with 
transmission amplitudes 2^1^/(1 + w\ 2 ) ( see 

is the modified spectral determinant. The matrix involved in the current density in the arc a 
is : 

t d£ _ 2 sin f 2w lW2 s b ( -£ 21 e^ E^e"*- \ , . 

from which we get the contribution of the scattering state ^^(x) to the current density in the 
arc a : 

= — ( E 1 "-^- | = —L- \-wf sin 6 sin kl + w? wl(st + s a s b cos 0)) . (83) 
2i7r V Ma/n tt\S\ 2 11 1 2V b n V ' 

Let us now study the weak coupling limit w\ 2 — ► 0. Close to a resonance, we obtain : 
where 7 = w ^ w ^ t Integrating the contribution of the resonance peak, we get : 



, dk2 kj^Hk?)*^^%Un-f) . (85) 
k±-8K wf + w^ V \ 2vr J 

In the r.h.s we recognize the persitent current of the level of the isolated ring — -^g(k^) 2 = 
~W ( 2?r ; ) ! multiplied by the "relative weight" w ™+ w -± of the scattering state ^ l \x), as 
expected from (|76|1. 

4.6 Graphs with localized states 

In this section we discuss the consequence of the possible existence of localized states in certain 
graphs. These states are not probed by scattering, consequently their contributions to the 
current is not given by the expressions derived above. 

For the sake of simplicity our discussion will be focused on the example of a ring in the regime 
of the integer quantum Hall effect, with one edge state. The potential hill at the middle is called 
an antidot (figure |HJ left). This example must be thought more as a toy model to understand 
the idea of localized states in graphs, than as a reallistic model to describe current distribution 
in a quantum Hall device where the effect of screening is important (the interested reader will 
find some discussion on the nature of edge currents and the role of screening in [451 146j ) . The 
system can be modeled by a ring with chiral scattering at the vertices (figure |H1 right). The ring 
has two bonds, i.e. four internal arcs : two arcs a and b carrying fluxes 9 a and 9b and the two 
reversed arcs denoted with a bar : a and b. The two leads are described by arcs 1 and 2. 



19 





Figure 8: Left : A mesoscopic device in the regime of the IQHE with an antidot at the middle. 
One edge state is open. Right : The graph that models this arrangement. The scattering at the 
vertices is chiral. 



In the basis of arcs {a, b, a, b\l, 2}, its vertex scattering matrix and bond scattering matrix 
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(86) 



The off diagonal blocks are not anymore transposed due to the breaking of the time reversal 
symmetry at the vertices (chiral scattering), however we keep the same notation as above for 
simplicity since there is no possible confusion. On the other hand 
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jtclb+i6b 
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(87) 



We recall that Qij is the transmission amplitude from arc j to arc i due to vertex scattering, 
and Rij describes bond scattering due to the potential. 



Scattering. The scattering matrix can be computed from Q and R with equation 
the result is obvious here, due to the absence of multiple scattering : 



however 







jkl b +i6 b 





The current density in the arc a is given by the matrix : 



i 




(89) 



From (|55|) we get the contribution of the scattering state ip^(x) to the current density in the 
arm a : ja = ^r, whereas the contribution of ijj^(x) obviously vanishes ja 
latter scattering state does not send current into the arc a. 



0, since this 



Localized states. We follow the discussion of [111 I13j : if localized states are present, their 
discrete spectrum is given by solving det(^ f - Q int ) = 0. Here, we see that the equation indeed 
possesses a set of solutions since det(R^ — Q mt ) = e~ lkl (e~ lM — e~ ld ) where 9 = 9 a + 0j. The 
spectrum of localized states is k n = (2Trn+0)/l for n G N if 6 G [0, 2tt[, since k ^ by convention. 
These states describe a clockwise motion of the electron in the loop of the graph (right part of 
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figure IBJ). In the quantum Hall ring picture, they correspond to states whose wave functions 
are localized on the edge of the antidot (left part of figure |H1) • The current associated to the 
state of energy k\ in the arc a is — J^k^ = —2(2mr + 0)/l 2 . Note that if one introduces some 
scattering on the bonds, the localized states are hybridized with the states of the continuum 
and the discrete part of the spectrum disappears. 

The discrete spectrum also brings some contribution to the current in the arms of the ring, 
which cannot be obtained from the scattering properties. Since the state ip^ does not contribute, 
the total average current in arm a finally reads : 

rco 00 9 

J a = dEh{E)j^){E)-Y,hnt{k 2 n ) ^(2nn + e) (90) 

where fi(E) is the Fermi distribution for the lead 1 and f[ n t(E) the Fermi distribution for the 
localized states inside the graph. 

Now we can give the general expression for the current in the arc a for a graph with localized 
states. Since the discrete spectrum of localized states {E n } is given by the equation 

= det^ - Q int ) oc ~[[(E-E n ) (91) 

n 

we have 

- E/W)i(*«f) ( 

+ fdEf int (E)-Im^-ln(det(E)-Q M ) V (92) 

where fi n t(E) is the Fermi distribution associated to localized states. The first term is the 
contribution of the scattering states whereas the second is the contribution of the localized 
states. 



5 Summary 

In this paper we have studied the two first cumulants of the charge of a graph connected to 
infinite wires, as well as the distribution of currents in the wires inside the graph. In particular, 
we have shown the relation with the scattering matrix, allowing to study these quantities in an 
out of equilibrium situation, when the graph is connected to wires put at different potentials. 
We have obtained a formula for the average current and the current correlations that generalizes 
previous results known for the equilibrium situation |17( I29j . 

We have also emphasized that the scattering matrix contains information only on the con- 
tinuous part of the spectrum related to scattering states. If some states remain localized in 
the graph, they give an additional contribution to the current not taken into account by the 
scattering approach. 

We have considered the case of graphs weakly coupled to the leads. It is interesting to 
remark that the results obtained in this context are expected to be of much more generality 
than graphs, since the starting point was to use an approximation of the scattering state near 
a resonant level (|21|). a form of great generality. In particular, the contribution of the resonant 
level n to the average of some quantity X defined inside the graph reads 

(X(t)>("> ~ X n (Q(t))(») , (93) 
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where X n = ((f n \X\ tp n ) is the expectation of X in the eigenstate \ <p n ) of the isolated system. 
Similarly the contribution of the n-th resonant level to the correlations of two observables X 
and Y reads : 

S$>r(u,)~X n Y n S${u). (94) 

These results apply to a situation with narrow resonances (r n <C \E n+ i — E n \). We repeat 
that we have not considered the effect of electron electron interactions in this article (weakly 
connected devices with resonant tunneling present in principle Coulomb blockade). It would be 
interesting to incorporate some effects of interaction. This could be already done in a mean field 
approximation to describe the effect of screening in the charge and current distribution following 
Biittiker's approach [31121137]. 
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A Relation between the different unsymmetrized correlators 

Let us consider A and B two hermitian operators associated to physical quantities in our system. 
There are two unsymmetrized correlation functions: 

/+oo 
dre iuJT ((A(t + r) B(t)) - (A(t + r)> (95) 
-oo 

/+oo 
dr e iwT «B(t) A(t + r)> - (B(t)) (A(t + r)» . (96) 
-oo 

These correlators do not depend on t for a system in a stationary state 5 . Note that, even if we 
consider out of equlibrium situations, only stationary states are considered in the present paper. 

In this case, using time translation invariance of one and two point correlation fonctions, we 
have : 

/+oo 
dr ((B(t - r) A{t)) - (B(t - r)> (A(t))) (97) 
-oo 

and this gives: 

SaA") = SbA- u ) ■ ( 98 ) 
When A = B, it shows that one unsymmetrized correlator determines the other one: 

S a ,aH = SaA-u)- ■ (99) 

Finally, note that the correlator Saa^) 1S re& l since the correlator in time obeys SaA t )* = 
SaA~ t )- 

B Structure of the stationary states near a resonance 

If we consider a graph weakly coupled to the leads, we expect the stationary scattering states 
to be closely related to the eigenstates of the isolated graph. The purpose of the appendix is 
to demonstrate the precise relation. The relations we will obtain are very reminiscent of the 

5 Glassy systems are an example for which this is not possible since the system never reaches a stationary state 
(weak ergodicity breaking). 
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Hamiltonian approach of chaotic scattering [37j (see also |48j). In this latter case some small 
couplings are introduced between an isolated system and the leads, whereas we rather start from 
a situation where the coupling can be arbitrary large and study the weak coupling limit to see 
how the properties of the isolated graph emerges from its scattering properties. 

Let us consider a graph G, whose spectrum is supposed to be non degenerate for simplicity 
(the occurence of degeneracies leads to complications related to possible existence of localized 
states in the graph non probed by scattering I13j). In a first step we describe how the 
eigenstates of the Schrodinger operator in the graph are constructed and in a second step we 
will establish the relation with stationary states in the weak coupling limit. 

Isolated graph. The spectrum of the Schrodinger operator in the graph is given by the equa- 
tion : 5(7) = 0, where £(7) = EL (7 + ^n) is the spectral determinant, whose construction 
is explained in |411 I42| for the case of free graphs, in |431 149| for graphs with potential and in 
|44j for graphs with general boundary conditions (more general than the continuity of the wave 
function at vertices). If the wave function is continuous at vertices : 

S{l)=l V/2 II f^(a)V 1 detM( 7 ) . (100) 

(a/3) K(XXa P J 

The product runs over all the bonds of the graph. We recall that the functions f a /3( x ) involved 
in M(j) are the solutions of the Schrodinger equation on the bond [7 — d^ + V^ a ^(x)]f(x) = 
for an energy E = —7. In general ^(7) = possesses the same set of solutions as 

detM( 7 )=0. (101) 

We do not discuss here the case where the sets of zeros of both equations do not coincide, which 
is a little bit pathological and would require to refine the following arguments. Let us however 
quote few examples of free graphs (V(x) = 0) for which it is the case : the graph made of one 
line (in this case detM = 1 is independent of 7), the complete graph |421 111 j ... . 
The component of the wave function (p n (x) on the bond (a/3) is : 

V>n(a/3)(x) = e [A ^ x (ip n , a f af3 (x) + <p n>p e~ w ^ f Pa (z)) (102) 

where (p ntCe is the wave function at the vertex a and A a p = 6 a p/l a p the vector potential. (Do 
not confuse the label n of the eigenstate with the greek labels that designate vertices). If we 
gather the wave function at the nodes in the ^/-dimensional column vector ip n , the eigenstate of 
energy E n is solution of 

M(-E n )ip n = . (103) 



Normalization. The normalization condition for the eigenstate reads : 



f dx \(f n (x)\ 2 = [ dx \<Pn(af3)(x)\ 2 = 1 

J Graph {a/3) J0 



(104) 



If we use the following relations [IHj : 



'a/3 



' dx al 3 



dx a/3 f a /3(x af 3) 2 = -dy^j^-(a) (105) 



a/3 



dx af3 f af3 (Xal3)ff3a(Xal3) = 5 7~[~(^) ( 106 ) 



a/3 
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we obtain 



/ dx |^(*)| 2 = £ U a 5 7 (-^(a)) <p n ,a +<Pn,ad 7 (^{(3)^ A <p n ,p 
J Graph L V dx a/3 / \ d:c a/3 / 



(107) 



If we replace the sum over bonds by a sum over vertices, the matrix M appears. Finally, the 
normalization condition reads for the V^-vector <p n : 

<piMy/7M(rf)]<p n = 1 , (108) 

where the spectral parameter is taken, after derivation, equal to the eigenenergy 7 = — E n — i0 + . 

Graph weakly connected to leads. When the graph is weakly coupled to leads (w a — ► 0) 

we expect that the stationary state ijjg (x) is proportional to the wave function of the isolated 

graph near the resonance E ~ E n : ift^ (x) oc <p n (x) for x £ Q. The question is how to recover 
precisely this relation from our formalism ? 

The resonance width. As a preliminary question, it is instructive to find an expression for the 
resonance widths. For this purpose, let us consider the determinant of the scattering matrix 

, „ , Nf det(M- W T W) , , 

and find an approximation near a resonance E n . 

For any fixed energy E > 0, M is an antihermitian matrix and can be written in terms of 
its purely imaginary eigenvalues \\ a (E) and its associate eigenvectors v a (E) : 

v 

M(-E) = iY J ^a(E)v a (E)vi(E). (110) 
0=1 

The eigenvectors are normalized as VaV a = 1. If the energy E is equal to the energy E n of 
an eigenstate of the isolated graph, one of the eigenvalues of M is vanishing : \i(E n ) = 0. 
We suppose the spectrum of the isolated graph to be non degenerate. The eigenvector v\(E n ) 
coincides with the eigenstate : v\{E n ) = u~ 1 ip n ; however these vectors are not normalized in 
the same way and differ in the multiplicative factor u n . 

Since det M(—E) is proportional to the spectral determinant and the spectrum supposed 
to be non degenerate, the eigenvalue \\(E) behaves linearly near the energy E n : X\(E) ~ 
(E — E n )/3 n . The normalization condition Q1U8|) reads : 

/ v \ 



t d E -iVEY,^a(E)v a (E)vUE) 



0=1 



, , _ , E=E„ 

then 



Vn = 1 (111) 



- P n ^/E~ipiv 1 (E n )v\(E n )i Pn = l (112) 

We obtain the normalization constant : v n = l/y/—k n (3 n where E n = k^. 

We now come back to detS. In the weak coupling limit w a —* we can compute perturba- 
tively the eigenvalues of M ± W T W to express the determinant : 

V 

det(M{-E) ± W T W) ~ \\(\\ a {E)±vi(E)W rv Wv a {E)\ (113) 



E~E,. 

a=l 



V 



ifin (E - E n ) ± vliE^W^Wv^En)) \{ i\ a (E n ) (114) 



a=2 
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We can use the relation v\{E n ) = y/—k n (5 n (p n to get : 

det(M(-k 2 )±W T W) ~ (k-k n ±^ n W T W Vn \2ip n k n l[iX a (k 2 n ). (115) 

Then : 

detSoc k ~l n ~ hn (116) 
k k n + vy n 

where the resonance width in fc-scale is : j n = ^ipnW T Wip n = | Y2 a =l w V\ i Pri,a\ 2 ■ This result 
is very satisfactory since it shows that the lead a brings a contribution to the resonance width 
proportional to the transmission probability w\ between the graph and the lead 6 and to the 
probability density |v?n,«| 2 = |(/? n (a)| 2 associated to the eigenstate of the isolated graph, taken 
at the vertex a where the graph is connected. In energy scale the resonance width reads : 



T n = 2k nln = ^/E~ n ^ n W T W^ n = T n,a (117) 

a=l 

where 

r„, Q = ^Kw 2 a \ip n {a)\ 2 (118) 



is the contribution of the lead a. 

The wave function. We recall that the V x L-matrix that gathers the values of the L scattering 
states at the V vertices is ^U] : 

^^swTwm wT - (119) 

We call the V-column vector gathering the values of ip^ a \x) at the V vertices : ip^ = 

■ • • , if)y^) T . The matrix \& is obtained by gathering these L column vectors : \l/ = 

■ ■ ■ , ip ^). In the weak coupling limit (w a — > 0) and near the resonance E n we can 
keep only the contribution of the vanishing eigenvalue \\i{E) of M to compute : 



-r vi(E n )v\(E n ) . (120) 



M + W^W if3n (E - E n ) + v\(E n )W^Wv 1 (E. l 
It follows that the scattering state at the vertex /x is : 

7(a) vr, 1 i/ 2 , tw/T\ 1 iw ^ { Pn,a 



/^fe A: - A:„ + i7„ ' n ^4nk k - k n + 

Since the vertex fi could be any point of the graph because we have always the freedom to 
introduce an additional vertex of weight A = on any bond without changing the properties of 
the graph, we can rewrite more elegantly : 



Or using energy scale : 



= *™fi a) W.W (122) 
k^k n V^k n k- k n + ry n 



£~£„ E - E n + iF n 



s the transmission amplitude is 1 2 ^ a for finite w a |10) 
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The local density of states. The contribution of the scattering state ip^ (x) to the off-diagonal 
LDoS (x \5{E-H)\ x') is : 

#(*)$?V) - (F (124) 
The LDoS is obtained by summing these contributions over a. 

T/ie scattering matrix. The same discussion can be done to find an approximation for the 
scattering matrix. We obtain the well-known Breit-Wigner structure : 

Y, a p(E) ~ -d al 3 + • (125) 

E~E n Hj — £j n + ll n 

C Matrix elements of the charge operator 

Our aim is to relate 

p^(E) = j dx^\xf ^\x) = £ dx^Ux) jW(x) (126) 

J Graph UjuA ^ 

to the scattering matrix. The sum runs over the B bonds. The relation (|14j) was proven for 
a = (3 in 111] using a different method. 

The computation of this integral follows exactly the lines of the one done discussing the 
normalization of the states in the isolated graph. Then we obtain : 



p W(E) = ffi* 9e (iy/EM^) ffl , (127) 



that is 



p {<*,$){E) = (§tA(i/EM)^) (128) 

a/3 



dE 



■Uw- 1 



2ivr V -M + W T W 
where we have used ()119|) . From © and 



dM 1 „, 

2 + — M 

dE E 



W T ) (129) 



M + W T W ; ap 



™ = _ 2W 1 dM 1 T 

dE M + W T W dE M + W T W y ' 



we finally obtain the desired relation : 



(131) 



An alternative way to relate p( a '^ (E) to derivative of the scattering matrix is to introduce 
the variable conjugate to the charge : a constant potential U in the graph. The total potential 
now reads V(x) + U Og(x) where 9g(x) = 1 if x € Q and 6g(x) = if x belongs to the leads. In the 
presence of U, the function f a p(x) involved in M is solution of [E + d^ — V( Q| g) (x) — U]f a p{x) = 0. 
These functions are obtained by a shift of the spectral parameter : f^Jx;E) = f®Jx;E — U). 
It immediately follows that : 



pt<x,P)( E ) = -[¥—( iVem) . (132) 

V dU V J J a/3 
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Using the same arguments as above we get 



^--55 (*§)„■ < 133 > 

It is interesting to compare this relation with (|131|) : it shows that it is not similar to differentiate 
with respect to a constant potential U or with respect to the energy since the potential does not 
live in the wires. The difference however vanishes at high energy (WKB limit). 

A relation between the stationary states and the functional derivative of the scattering matrix 

was proven for graphs ^3] where it is explained how it can be computed with algebraic calcu- 
lations only. It follows that we can rewrite the equation (|131|) : 

«5£ dS dS 1 , 7 x , . 

The first equality, which is obtained by identification of (j!26|134j) with (|133j) . is also a conse- 
quence of the definition of the functional derivative. 
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